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LEBESGUE CONSTANTS ARISING IN 
A CLASS OF COLLOCATION METHODS * 

WILLIAM W. HAGERt, HONGYAN HOUt, AND ANIL V. RAO§ 


Abstract. Estimates are obtained for the Lebesgue constants associated with the Gauss quadra¬ 
ture points on (—1, -i-1) augmented by the point —1 and with the Radau quadrature points on either 
(—1, H-l] or [—1, +1). It is shown that the Lebesgue constants are 0(v7V), where N is the number of 
quadrature points. These point sets arise in the estimation of the residual associated with recently 
developed orthogonal collocation schemes for optimal control problems. For problems with smooth 
solutions, the estimates for the Lebesgue constants can imply an exponential decay of the residual 
in the collocated problem as a function of the number of quadrature points. 
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1. Introduction. Recently, in [Sj ID 0 El HOl HU EO], a class of methods was 
developed for solving optimal control problems using collocation at either Gauss or 
Radau quadrature points. In m and m an exponential convergence rate is es¬ 
tablished for these schemes. The analysis is based on a bound for the inverse of a 
linearized operator associated with the discretized problem, and an estimate for the 
residual one gets when substituting the solution to the continuous problem into the 
discretized problem. This paper focuses on the estimation of the residual. We show 
that the residual in the sup-norm is bounded by the sup-norm distance between the 
derivative of the solution to the continuous problem and the derivative of the inter- 
polant of the solution. By Markov’s inequality [18], this distance can be bounded in 
terms of the Lebesgue constant for the point set and the error in best polynomial 
approximation. A classic result of Jackson m gives an estimate for the error in best 
approximation. The Lebesgue constant that we need to analyze corresponds to the 
roots of a Jacobi polynomial on (—1,-|-1) augmented by either r = -|-lorr = —1. 
The effects of the added endpoints were analyzed by Vertesi in [24]. For either the 
Gauss quadrature points on (—1, -|-1) augmented by r = -1-1 or the Radau quadrature 
points on (—1, -|-1] or on [—1, -|-1), the bound given in [24] Thm. 2.1] for the Lebesgue 
constants is 0{\og{N)'/N), where N is the number of quadrature points. We sharpen 
this bound to 0{'/N). 

To motivate the relevance of the Lebesgue constant to collocation methods, let 
us consider the scalar first-order differential equation 

x{t) = f {x{t)) , rG[-l,-bl], a;(-l) = xo, (1.1) 

where / : R —>■ R. In a collocation scheme for (HI]), the solution x to the differential 
equation HU is approximated by a polynomial x that is required to satisfy the 
differential equation at the collocation points. Let us consider a scheme based on 
collocation at the Gauss quadrature points —1 < ti <T 2 < ... < <+1, the roots 
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of the Legendre polynomial of degree N. In addition, we introduce the noncollocated 
point tq = —1. The discretized problem is to find x £ Vn, the space of polynomials 
of degree at most N, such that 

xijk) = f{x{Tk)), l<k<N, a;(-l) = a;o. (1.2) 

A polynomial of degree at most N is uniquely specified by TV + 1 parameters such 
as its coefficients. The N collocation equations and the boundary condition in (II.2p 
yield TV + 1 equations for the polynomial. 

The convergence of a solution of the collocated problem (O to a solution of the 
continuous problem (El) ultimately depends on how accurately a polynomial inter- 
polant of a continuous solution satisfies the discrete equations (US). The Lagrange 
interpolation polynomials for the point set {tq, ti, ..., tn} are defined by 

N 

L^{T) = Y[^^, 0<i<N. (1.3) 

j=o 

The interpolant x^ of a solution x to (ED is given by 

N 

X^{t) = ^x{Tj)Lj{T). 

3=0 

The residual in (11.21) associated with a solution of (11.11) is the vector with components 

ro = x^(-l) - lo, Tk = x^{Tk) - f{x^{Tk)), 1 < fc < TV. (1.4) 

For the Gauss scheme, rg = 0 since x satisfies the boundary condition in (11.11) . The 
potentially nonzero components of the residual are rt, 1 < fc < TV. 

As we show in Section [2l the residual can be bounded in terms of a Lebesgue 
constant and the error in best approximation for x and its derivative. The Lebesgue 
constant Ajv relative to the point set {tq, ti, ..., tat} is defined by 

Aat = max|y]|Lj(T)| : r G [-1,1]| . (1.5) 

The article [1] of Brutman gives a comprehensive survey on the analysis of Lebesgue 
constants, while the book [19] of Mastroianni and Milovanovic covers more recent 
results. 

The paper is organized as follows. In Section [S] we show how the Lebesgue con¬ 
stant enters into the residual associated with the discretized problem (11.21) . Section [3] 
summarizes results of Szegd used in the analysis. Section |4| analyzes the Lebesgue 
constant for the Gauss quadrature points augmented by r = — 1, while Section [5] 
analyzes Radau quadrature points. Finally, Section [51 examines the tightness of the 
estimates for the Lebesgue constants. 

Notation. Vn denotes the space of polynomials of degree at most TV and || • || 
denotes the sup-norm on the interval [—1,-|-1]. The Jacobi polynomial 
TV > 1, is an TV-th degree polynomial, and for fixed a > — 1 and /3 > — 1, the 
polynomials are orthogonal on the interval [—1,-|-1] relative to the weight function 
(1 — t)“( 1 -|- r)^. Pat stands for the Jacobi polynomial or equivalently, the 

Legendre polynomial of degree TV. 
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2. Analysis of the residual. As discussed in the introduction, a key step in the 
convergence analysis of collocation schemes is the estimation of the residual defined in 
(na). The convergence of a discrete solution to the solution of the continuous problem 
ultimately depends on how quickly the residual approaches 0 as A tends to infinity; for 
example, see Theorem 3.1 in [5, Proposition 5.1 in [12], or Theorem 2.1 in m- Since a 
solution X of (El) satisfies the differential equation on the interval [—1, +1], it follows 
that x{Tk) = f{x{Tk)), 1 < k < N. Hence, the potentially nonzero components of the 
residual can be expressed — x{Tk), 1 < k < N. In other words, the size 

of the residual depends on the difference between the derivative of the interpolating 
polynomial at the collocation points and the derivative of the continuous solution 
at the collocation points. Hence, let us consider the general problem of estimating 
the difference between the derivative of an interpolating polynomial on the point set 
tq < Ti < ... < Tvr contained in [—1, +1] and the derivative of the original function. 

Proposition 2.1. If x is continuously differentiable on [—1,+1], then 

lli — x^ll < (l + 2A^) inf lli; — g|| 

+N'^{1 + An) inf ||a:-_p|| (2.1) 


where x^ G Vn satisfies x^{xk) = x{Tk), 0 < k < N, and Ajy is the Lebesgue constant 
relative to the point set {tq, ti, ..., tn}- 

Proof. Given p G Vn, the triangle inequality gives 

||i; — < ||i: — pII + ||p — x‘^11 . (2.2) 

By Markov’s inequality [18], we have 


Ip — x^\ 


< 

< N^An max \p(Ti) 
0<i<N 


N 

^(p(Ti) -x(tO)L*(t) 


2=0 

xiTi)\ < A^AtvIIp- a:||. 


(2.3) 


Let q G Vn with (/(—I) = a;(—1). Again, by the triangle and Markov inequalities, we 
have 


lli-pll < Hi - g|| + ||9-p|| < ||i- g|| +N'^\\q-p\\ 

< l|i-9ll+^^(lk-a;|| + ||a;-p||). (2.4) 


By the fundamental theorem of calculus. 


\q{t) - X 



(?(«) 


i(s)) ds 



i(s)| ds < 2\\q — i||. 


We combine this with (12.41) to obtain 

Ili-Pll < (1 + 2A2 )||x - q|| +N'^\\x-p\\. (2.5) 

To complete the proof, combine (12.2L (12.31) . and (12.51) and exploit the fact that 
{q ■■ ^(-l) = a;(-l), q G Vn} = {q : q G Vn] ■ 


□ 
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An estimate for the right side of (12.11) follows from results on best uniform ap¬ 
proximation by polynomials, which originate from work of Jackson m- For example, 
the following result employs an estimate from Rivlin’s book pT| . 

Lemma 2.2. If x has m derivatives on [—and N > m, then 


inf 

P&Vn 


Ik-Pll < 


12 



where denotes the m-th derivative of x. 
Proof It is shown in [21] Thm. 1.5] that 


inf lla; 

P&Vn 


pII < 


6 


m -I- 1 



1 


N — m 


( 2 . 6 ) 


(2.7) 


where iVm is the modulus of continuity of x^'^\ By the definition of the modulus of 
continuity, we have 


1 


N — m 


= sup 




■■ Ti,T2 G [-1,1], |ti 


T 2 I < 



Since 


|x('")(ti)-x('")(t 2)| <2||a;W|l, 


(12.61) follows from (12.71) . □ 

If Aat = 0{N) and m > 4, then Proposition 12.11 and Lemma 12.21 imply that the 
components of the residual approach zero as N tends to infinity. Moreover, if x is 
infinitely differentiable and there exists a constant c such that ||a::^’"^|] < c™, then we 
take m = N — 1 m Lemma 12.21 to obtain 


inf 

peVN 


Ik 



Hence, the convergence is extremely fast due to the 1/7V^ factor. 

3. Some results of Szego. We now summarize several results developed by 
Szegd in [22] for Jacobi polynomials that are used in the analysis. The page and 
equation numbers that follow refer to the 2003 edition of Szegd’s book published by 
the American Mathematical Society. First, at the bottom of page 338, Szegd makes 
the following observation: 

Theorem 3.1. The Lebesgue constant for the roots of the Jaeobi polynomial 
P^’^\t) is 0(iV°-®+''') if j := max(Q;,/3) > —1/2, while it is 0{logN) if ^ < —1/2. 


For the Gauss quadrature points, a = /3 = 0, 7 = 0, and kjq = 0{\fN). The 
result that we state as Theorem 13.11 is based on a number of additional properties of 
Jacobi polynomials which are useful in our analysis. The following identity is a direct 
consequence of the Rodrigues formula [^ p. 67] for P^’^\ 

Proposition 3.2. For any a and /3 G R, we have 


’“^(-'t) for all T G [-1,+!]. 


(3.1) 
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The following proposition provides some bounds for Jacobi polynomials. 
Proposition 3.3. For any a and /3 G K and any fixed constant ci > 0, we have 


P^“’^^(cos0) = < 


OiN°‘) z/6ie[ 0 

6»-“-o-5o (iv-i/2) ifee[ ciiv-i 

(tt-(A f-1/2) z/6»e[ tt /2 

0{N^) z/6ie[ TT-ciiV-i 


ciN-^ 
■K12 

TT — 

TT 


Proof. The bounds for 9 G [0,cA^ and for 6 € [cN ^, 7 r/ 2 ] appear in 
(7.32.5)]. If 0 G [ 71 / 2 , tt], then tv — 6 G [0,7r/2] and by (13.11) . 


P^“’^^(cos 6 ») = cos( 7 r - 6 »)) = (-l)^P^^’“^(cos( 7 r - 6 »)). (3.2) 


Hence, for 6 G [ 7 r/ 2 , 7 r], the first two estimates in the proposition applied to the right 
side of (13.21) yield the last two estimates. □ 

The next proposition provides an estimate for the derivative of a Jacobi polyno¬ 
mial at a zero. 

Proposition 3.4. If a > —1 and fi > —1, then there exist constants 72 > 7i > 0, 
depending only on a and fi, such that 


71* 


-P-l.bj^P+2 ^ 




(.Oi.P), 


N 


in) < 


whenever < 0 where ri < r 2 < ... < tjv are the zeros of (the smallest zero 

is indexed first). Moreover, if 9i G [ 0 , 7 r] is defined by cos9i = Ti, then there exist 
constants 74 > 73 > 0 , depending only on a and fi, such that 




p(“./3) 



< JiVNiTT - 9^)-<^-^■^ 


(3.3) 


whenever 9i G [ 7 r/ 2 , 7 r]. 

Proof. In (8.9.2)], it is shown that there exist 72 > 7 i > 0, depending only 
on a and /3, such that 


7 ^ .-/ 3 - 1 . 5 ^/ 3+2 < 


P 


(/3.a) 


N 


(C7z) 


< 72* 


-/ 3 - 1.5 


JV^+2 


(3.4) 


whenever Ci > 0 where cti > tT 2 > • ■ • > ctn are the zeros of P^’“^ (the largest zero 
is indexed first). By Proposition 13.21 Ti is a zero of if and only if —Ti is a zero 

of P^’“^ Hence, the zeros of are —ri > —T 2 > ... > —tn- Moreover, 

P^“’'')(r) = ±Plf’“)(-T). (3.5) 

The bound given in the proposition for \Pj^’^\Ti)\ with Tj < 0 is exactly the bound 
(13.41) for |P^^’“^(cri)| with ai > 0. 

It is shown in (8.9.7)], that there exist constants 74 > 73 > 0, depending only 
on a and /3, such that 









(3.6) 


whenever fii G [ 0 , 7 r/ 2 ] where cosfii = ai. Since cosfii = ai = —Ti = cos( 7 r — 0i), it 
follows that fii = TT — 9i, and (13.51) and (13.61) yield (13.31) . □ 
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4. Lebesgue constant for Gauss quadrature points augmented by — 1. 

In this section we estimate the Lebesgue constant for the Gauss quadrature points 
augmented by tq = — 1. Due to the symmetry of the Gauss quadrature points, the 
same estimate holds when the Gauss quadrature points are augmented by +1 instead 
of —1. The Gauss quadrature points are the zeros of the Jacobi polynomial 
which is abbreviated as By Theorem l3.ll the Lebesgue constant for the Gauss 

quadrature points themselves is 0{'/N)- The effect of adding the point rg = —1 to 
the Gauss quadrature points is not immediately clear due to the new factor (1 + r^) in 
the denominator of the Lagrange polynomials; this factor can approach 0 since roots 
of Pn approach —1 as TV tends to infinity. Nonetheless, with a careful grouping of 
terms, Szego’s bound in Theorem [ST] for the Gauss quadrature points can be extended 
to handle the new point tq = —1. 

Theorem 4.1. The Lebesgue constant for the point set consisting of the Gauss 
quadrature points —1 < ti < r 2 < ... < r^v < +1 {the zeros of Pm) augmented with 
Tq = — 1 is 0{\fN). 

Proof. Define 


l)r) = {t - ti){t - T2) ■. .{t - tn), and L(r) = (r + l)l(r). 
The derivative of L{t) at Ti is 

r ^(-1), * = o, 


L{Ti) = l{n) + in + i)l{n) = 


in + l)^('ri), i > 0. 


Hence, the Lagrange polynomials Liir) associated with the point set {tq, ti, ..., tn} 
can be expressed as 


Lz(t) = 


L{t) 


L(Ti)(T - n) 


/(r)A(-l), 

L(t) 


in + i)l{n)iT - n) 

Since Pm is a multiple of I (it has the same zeros), it follows that 

Pm{t)IPm{-1), i=0, 

(r + l)PAr(r) 


1 = 0 , 

, i > 0. 


(4.1) 


L,{t) = 


in + l)TW(Ti)(r - n) 


, z > 0. 


By gU (7.21.1)], |PAr(T)| < 1 for all r € [-1,+1], and by [HJ (4.1.4)], lPAr(-l)l = 
(—1)^. We conclude that \Lq{t)\ < 1 for all r G [—1,+!]. Hence, the proof is 
complete if 


N 


max |Lz(t) 1 : r G [-1,1] > = 0{Vn). 


(4.2) 


. i=l 


For any r G [—1, +1], the integers i G [1, N] are partitioned into the four disjoint 


11 = {i e [1, A^] : n > 0}, 

1 2 = {i & [1, iV] : -1 < Ti < 0, Ti > r}, 

1 3 = {i € [1, A^] : -1 < Ti < 0, Ti < T, T - Ti < Ti + 1}, 

1 4 = {* G [1, A^] : -1 < Ti < 0, Ti < T, T - Ti > Ti + 1}. 


sets 
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Let Ii 23 denote Ii U I 2 U I 3 . Observe that for any i G I 123 and r £ [—1,+1], 
(r + l)/(ri + 1) < 2. Consequently, for all i G X 123 , 


\L^iT)\ = 


(r + I)Pjv(t) 


in + l)PAr(ri)(r - Ti) 


< 


2|^^n(t)| 


\PN(.Ti){T -Ti)\ 


This bound together with Theorem 13.11 imply that 


2|Pv(rL < -= 0{VN) 


Y \Ht)\ <y _ -. - 

I^N(r.)(r - r,)| l-PN(r.)(T - n)\ 

since the terms in the final sum are the Lagrange polynomials for the Gauss quadrature 
points. To complete the proof, we need to analyze the terms in (14.2|) associated 
with the indices in I 4 . These terms are more difficult to analyze since + 1 in the 
denominator of Li could approach 0 while r + 1 in the numerator remains bounded 
away from 0 . 

For i G I 4 , we have 

T + 1 = (t - Ti) + (tj + 1) < 2(r - Ti) 
since t — ti > + 1. Hence, 

2 |Pw(t)| . 2 


\L^{T)\ < 


< 


|(Ti + l)P/v(rO| |(ri + l)PAr(Ti)| 
since |P/v('7')| < 1 for all r G [—1,+1] by [HI (7.21.1)]. It follows that 

Y\L^{T)\<Y ——■—< y" ——■—■ 

Given 0 G [ 7 r/ 2 , 7 r], define (p = it — 6. Observe that 


(4.3) 


1 + COS 0 


2cos^(6(/2) sin^((^/2) xe[o. 7 r/ 4 ] sin^ a; 4 


Hence, for 0 G [ 7 r/ 2 , 7 r], we have 


1 + cos 0 > 


= Y^-Y- 


(4.4) 


By the bounds [22l (6.21.5)] for the roots of the Jacobi polynomial when a and 

j3 G [—0.5,+0.5], it follows that 


2i- 1 


TT < TT — 0,- < 


2i 


2N+1 


2N + 1 

where cos0i = Ti. This implies the lower bound 

2i-l 


TT -0i> 


2N + \ 

We combine (14.41) and (14.61) to obtain 


TT > 


i 

m 


TT, 1 < * < 


N' 


1 + T. > -^{7T-0if > -y Ur 


(4.5) 


(4.6) 


(4.7) 
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By Proposition inm 

|P/v(cos6»i)| > 

This lower bound for the derivative and the lower bound (14.71) for the root imply that 

1 


< 




,-1/2 


{1 + Ti)\pN{Ti)\ \ 471 , 

Hence, we obtain the following bound for the I 4 sum in (14.31) : 


E 

— l<ri<0 


2 

(1 + Ti)\PN{Ti)\ 




This bound inserted in (j4.3l) completes the proof. □ 

5. Lebesgue constants for the Radau quadrature points. Next, we esti¬ 
mate the Lebesgue constant for the Radau quadrature scheme. There are two versions 
of the Radau quadrature points depending on whether ti = — 1 or tjv = +1- Since 
these two schemes have quadrature points that are the negatives of one another, the 
Lebesgue constants are the same. The analysis is carried out for the case tn = -|-1. 
In this case, the Radau quadrature points are the TV — 1 roots of augmented by 

Tjv = 1. Szego shows that the Lebesgue constant for the roots of is 0(TV^/^). 

We show that when the quadrature point t/v = 1 is included, the Lebesgue constant 
drops to 0{VN). 

The analysis requires an estimate for the location of the zeros of Our 

estimate is based on some relatively recent results on interlacing properties for the 
zeros of Jacobi polynomials obtained by Driver, Jordaan, and Mbuyi in [^. Let r^' 
and r", z > 1, be zeros of Pn-i and P/v respectively, arranged in increasing order. 
Applying [51 Thm. 2.2], we have 


<Ti < T,- 


z = 1, 2,..., TV — 1, where —1 < ri < T 2 < ... < tn-i < +1 are the zeros of • 
Let 9i G [0,7r] be dehned by cos^i = n. By the estimate (14.51) for the zeros of Pat, it 
follows that the zeros of Pj^fl have the property that 


2 t-l\ ^ /2(T + 1)' 

TT < -- I tt, 


2TV - ly " ' ‘ ' \ 2N+1 

When i is replaced by N — i, these bounds become 


2 i- 1 


TT < TT — 0i < 


2 i 


2N-1 


1 < z < TV - 1. 


1 < z < TV - 1. 


,2TV-b 

Together, (15.11) and (15.21) imply that 

TT — 6 i > i/N and > i/N, 1 < z < TV — 1; 


(5.1) 


(5.2) 


(5.3) 


moreover, taking into account both the upper and lower bounds, we have 

'lOTV- 7" 


, , /4(z + TV) + 2TV+l ^ 

+ V 4A12 _ 1 ' - 

f5(2N-l)\ 2.57r 

^ V 4A12 _ 1 


4A12 _ 

1 < z < TV - 2. 


(5.4) 
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Thus, the interlacing properties for the zeros leads to explicit bounds for the separation 
of the zeros; for comparison, Theorem 8.9.1 in yields 9i — 0^+1 = 0{1)/N, while 
(EH) yields an explicit constant 2.57r. These estimates for the zeros of are used 

to derive the following result. 

Theorem 5.1. The Lebesgue constant for the Radau quadrature points 


— 1 < Ti < r2 < ... < Tat = 1 

{the zeros of augmented by tjv = +1) is Ofs/N). 

Proof. The Lagrange interpolating polynomials Ri, 1 < i < N, associated with 
the Radau quadrature points are given by 


R^iT) = 


1 — T 


.1 - Ti 

Similar to (El), the Ri can be expressed 


1 1 

N-l 

l<i<N-l, Rn{t) = 

■ 1 ' * '7 

J = 1 

i=i 


- Ti 


- Ti)PN-l{n){T - n)' 


Ri{t) = < 


i < iV, 


i = N. 


(5.5) 


By E (4.1.1)] and [22l (7.32.2)], we have 


^jvl°i(l) = N and |T’j)/l\^(t) 1 < iV for all t G [-1,+1]. 


5 ( 1 . 0 ) 


(5.6) 


We conclude that ji?Ar(r)j < 1 for all r G [—1, +1]. Hence, the proof is complete if 


'N-l 


max < ^ |7?i(r)l : r G [-1, +1] > = 0{Vn). 


(5.7) 


i=l 


Let ^ > 0 be a small constant. Technically, any S satisfying 0<i5<l/2is small 
enough for the analysis. Szegd establishes the following bounds when analyzing the 
Lebesgue constants associated with the roots of Jacobi polynomials: 


N 


E 


p. 


( 1 . 0 ) 


N 


(t) 


-Pjv’°^(H)(r-ri) 


0{VN) if tG [-1,J-1], 
OflogN) if r G [J — 1,1 — J], 
0(iV3/2) if T-g [i_ j^i]. 


(5.8) 


Szegd considers the general Jacobi polynomials on pages 336-338 of [22], while 

here we only state the results corresponding to a = 1 and /3 = 0. 

We first show that (15.71) holds when r G [—1,1 —J]. Observe that (1 —t)/( 1 —Ti) < 
4/(5 when Ti < 1 — 5/2 and t G [—1, +1]. It follows from (15.81) that 




Ti<l-S/2 


E 

Ti<l-S/2 




p 


N-liji){T-‘ 


I 0{y/N). tG[-1,,5-1], 

\ O(logiV), tG [,5-l,l-5]. 


(5.9) 
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When Tj > 1 — (5/2 and r G [—1, +1 — (5], we have \t — Ti\ > (5/2; hence, 


E s (i 

i>n>i-<5/2 ' 


E 

l>ri>l~5/2 




{n - l)P^j^'°l{Ti) 


(5.10) 


We have the following bounds for the factors on the right side of (15.101) : 

(a) By Proposition 13.31 = 0{1) if r € [—1,5— 1] and = 

0{N-^/^) if tG [5-1,1-5]. 

(b) By (l?^ . - 1, where cosS* = > 0. 

(c) By a Taylor expansion around 9 = 0, 

0^/4:<l-cose <9^/2, 6»G[0,7r/2]- (5.11) 

By (b) and the lower bound in (c) at 0 = 0^, we have 

(1 - T.)|PE°i(r,)| > O.25730-^/"yiV^. (5.12) 


We combine this with (a) and (15.101) to obtain 


E =' 

l>ri>l-6/2 


N 


0(5V-i/2)^ 

N 

o(iV-i)^v^ 




= 0{Vn), tG [-1,5-1], 

= 0(1), tG [5-1,1-5], 


since 9i G [0,7r]. This establishes (15.91) for all r G [—1,1 — 5]. 

To complete the proof of (EH), we need to consider r G (1 — 5,1]. The analysis 
becomes more complex since Szego’s estimate (15.81) is 0{N^/‘^) in this region, while 
we are trying to establish a much smaller bound in (15.71) : in fact, the bound in this 
region is O(logiV) as we will show. For the numerator of Riir) and r G (1 — 5,1], 
Proposition 13.31 and (15.111) yield 


(1 - 'r)l-Fi-°i('r)| = (1 - cos6»)|P^^i°i^(cos6»)| 


9^0{N), 9€[0,1/N], 

6)1/2(9(7v-i/ 2)^ 6(G [l/5V,7r/2], 

(5.13) 


Given r G (1 — 5,1], let us first focus on those i in (15.71) for which jr — r^j > 5. In this 
case, Ti < 1 — 5 or 1 — Ti > 5, and (15.131) gives 


E 

\T — Ti I >6 


= E 


|t —Ti I >(5 


(r - 1 )pE°i^(t) 


{n - - Ti) 




52 


\T_T^\>S i-PEi(p)I 


(5.14) 


The lower bounds (13.31) and (13.61) imply that 

^ |P.(r)|=O(iV-i)^0,®/" + O(5V-i)^ |7r-0,|3/2 = O(l), (5.15) 

\T—Ti\>s n>o Ti<o 
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since the terms in the sums are uniformly bounded and there are at most N terms. 

The next step in the proof of (15.71) for r £ (1 — 1] is to consider those terms 

corresponding to |t —Ti| < 6. Since i5 is small, it follows that both r and are near 1, 
and consequently, 0 and 0i are small and nonnegative, where cos 9 = r and cos 6i = r^. 
In particular, 0 < 9i < t: 12. In this case where Ti is near t, it is important to take 
into account the fact that t — Ti is a, divisor of the numerator (t) . To begin, we 
combine the lower bound in (13.61) and the bounds in (15.111) to obtain 


(l-r) 


< 


202 






(5.16) 


It follows from (|5.5D that 








(5.17) 


The mean value theorem and the formula [HI (4.21.7)] for the derivative of 
P^-i{t) in terms of gives the identity 


(5.18) 


pU.U) 


pSl%) - pU.V.) 

fN + l\ 

P''^'ll{c0S'qi) 

T -Ti 


T -Ti 

V 2 ) 


where rji lies between 9 and 9i. Together, (15.171) and (15.181) imply that 

|i?,(r)| =O(iVi/2)02y^ P^Y2 ’(cos77,) . 


(5.19) 


The estimate (15.191) is useful when is near r. When Ti is not near r, we proceed 
as follows. Use the identity 

O o ■ (“ + /3) . (a-/3) 
cos a — cosp = —2 sm-^-sm-^-, 

to deduce that 


\t — Ti\ = \ COS9 — cos0j| > -tt\ 9^ — 9i 

'TT^ I 


(5.20) 


when \9 + 9i\ < tt, which is satisfied since both 9 and di are near 0. Exploiting this 
inequality in (15.171) yields 


|i?,(T)| =O(7V-l/2)02y^ 




02 - 02 


(5.21) 


Recall, that we now need to analyze the interval r £ [1 — d, 1] and those i for 
which |t — Ti\ < S. Our analysis works with the variable 0 £ [0,7r/2], where cos0 = t. 
The interval 0 £ [0,7r/2] corresponds to r £ [0,1] which covers the target interval 
[1 — S, 1] when S is small. By [221 (7.32.2)], we have 

\pI^L%0S7j,)\ < N{N - l)/2. 

If 0 £ [0, c/7V], where c is a fixed constant independent of N, then it follows from 
(15.191) that |i?i(T)| = Moreover, if 0^ < 20 < 2c/N, then |i?i(T)| = 0(1). 
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By the bounds (15.3p . the number of roots that satisfy < 2c/N is at most 2c, 
independent of N. On the other hand, if 6i > 26, then 9 < 9i/2 and 

| 02 _ 02 | ^Q2_g2 > ( 3 / 4 ) 02 ^ 

With this substitution in (j5.21[) . we have 


mr)\ 




p(l.O) 



By (15.6L (t)| < N. Hence, if 0 G [0,c/A^], then by (15.31) . we have 


N-l 

-3/2 _ nt^r-3/2^ Y- fl-3/2 


^ |i?.(r)| = 0(/V-3/2) ^ 0-3/2 = 0(/V-3/2) ^ 0 

\r — Ti\<5 \T — ri\<5 i 


N-l 


N-l 


= 0{N-^/^) (^N-! = 0(1) E = 0(1), 


for all 0 G [0, c/N]. 

Finally, suppose that 0 G [c/iV, 7 r/ 2 ]. By (15.41) the separation between adjacent 
zeros Oi and 0i+i is at most 2.5tt/N. Hence, if 6i is within k zeros of 0, then 77 ^ > 
0 — 7 = 2 . 57 rfe. Here fc > 2 is an arbitrary fixed integer. By Proposition 13.31 

we have 


Pn:!1 (cos rj,) 


(0-7iV-3)-O2o(7V-O2). 


Choose c > 27 . If 0 G [c/N,'k/2], then 0/2 > c/{2N) > j/N. Hence, 0 — 7 //V > 0/2 
and 


Pn:!1 (cos rji) 


(0/2)-3/2o(Ar-i/2) = 0-5/2o(jv-O2). 


Combine this with (I5.19P to obtain 

|i?,(r)| = 0 (l)vW^. 

when 0 G [c/7V, 7 r/ 2 ] and 0^ is within k zeros of 0. If 0^ < 0, then i?i(T) = 0(1). If 
9i> 9 and 9i is within k zeros of 0, then 0^ — 0 < 7 /iV, and 


di/9 < (0 + ')/N)/9 < 1 + 7 /c 


when 0 G [cA^, 7 r/ 2 ]. Thus |i?i(T)| = 0(1) when 0 G [cA^, 7 r/ 2 ] and 9i is within k zeros 
of 0 . 

This analysis of Ri when 0^ is close to 0 needs to be complemented with an 
analysis of Ri when 0^ is not close to 0 and 0 G [c/N,7r/2]. For 0 in this interval, 
Proposition 13.31 yields |P)^^G°^(cos0)| = By (15.211) . we have 


mr)\ 


0{N-^) 


V9y% 

|02 - 02 |' 


If 0 > 29i, then 02 - 02 > (3/4)02 


|i?,(T)| =O(Af-l)0-3/20l/2. 


(5.22) 
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By (15.2L we have 

|i?^_,(T)| =o((Af0)-3/2)V^TT. 


Recall that we are focusing on those i for which 0]v-i < 0/2. The lower bound 
0N-i > from (15.31) implies that i < N9/2 whenever 9N-i < 9/2. Hence, the set 
of i satisfying i < N9 is a superset of the i that we need to consider, and we have 


^ |i?,(r)|= ^ |i?^_,(r)| = O((iV0)-3/2) ^ Vm 

BiKSjl eiv-i<e/2 i<NB 

= 0{{N9)-^/'^){N9 + 1 ) 3/2 ^ 

On the other hand, if 0 < 20^, then we have 

102-0/1 “ i(0-0.)(0+0*)i - i(0-0z)0*i - i0-0.r 

Combine this with (I5.22p to obtain 


mr)\ 


0 ( 1 ) 

iiV0-iV0,r 


Earlier we showed that |i?i(r)| = 0(1) for those i where the associated 0i is within k 
zeros of 0. When 0^ is more than k zeros away from 0, we exploit the estimate (EB 
for the zeros to deduce that |A^0 —7V0i| behaves like an arithmetic sequence of natural 
numbers. Hence, the sum of the |i?i(T)| over these natural numbers, where we avoid 
the singularity, is bounded by a multiple of log N. This completes the proof. □ 

6. Tightness of estimates. At the bottom of page 110 in [24], Vertesi states 
some lower bounds for the Lebesgue function. In the case of the Gauss quadrature 
points augmented by tat+i = +1 and the Radau quadrature points with tn = +1, 
the associated Lebesgue function is of order '/N at t = (n + r2)/2, the midpoint 
between the two smallest quadrature points. It follows that the 0{\fN) estimates 
for the Lebesgue constant are tight. To study the tightness of the estimates, the 
Lebesgue constants were evaluated numerically and fit by curves of the form a\/l^ +6, 
10 < IV < 100 (see Figures 15/1116.2|) . A fast and accurate method for evaluating the 
Gauss quadrature points, which could be extended to the Radau quadrature points, 
is given by Hale and Townsend in m- Figure 16.1116.21 indicate that a curve of the 
form a-\/]V + 6 is a good fit to the Lebesgue constant. Another Lebesgue constant 
which enters into the analysis of the Radau collocation schemes studied in m is 
the Lebesgue constant for the Radau quadrature points on (—1,+1] augmented by 
To = —1. As given by Vertesi in [24l Thm. 2.1], the Lebesgue constant is O(logn). 
Trefethen [23] points out that the Lebesgue constant on any point set has the lower 
bound 


Aat > 


log(V)+ 0.52125..., 


due to Erdos [7] and Brutman [2]. For comparison. Figure l6/3l plots this lower bound 
along with the computed Lebesgue constant. When the number of interpolation points 
range between 10 and 100, the Lebesgue constant for the Radau quadrature points 
augmented by the point —1 differs from the smallest possible Lebesgue constant by 
between 0.70 and 0.84. 
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Fig. 6.1. Least squares approximation to the Lebesgue constant for the point set corresponding 
to the Gauss quadrature points augmented by —1 using curves of the form ay/N + b 



Fig. 6.2. Least squares approximation to the Lebesgue constant for the point set corresponding 
to the Radau quadrature points using curves of the form ay/N + b 



Fig. 6.3. Least squares approximation to the Lebesgue constant for the point set corresponding 
to the Radau quadrature points on (—1, +1] augmented by —1 using curves of the form a log N b 
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7. Conclusions. In Gauss and Radau collocation methods for unconstrained 
control problems pUfTS] . the error in the solution to the discrete problem is bounded 
by the residual for the solution to the continuous problem inserted in the discrete 
equations. In Section [2 we observe that the residual in the sup-norm is bounded by 
the distance between the derivative of the continuous solution interpolant and the 
derivative of the continuous solution. Proposition 12.11 bounds this distance in terms 
of the error in best approximation and the Lebesgue constant for the point set. We 
show that the Lebesgue constant for the point sets associated with the Gauss and 
Radau collocation methods is 0{'/N), and by the plots of Section [51 the Lebesgue 
constants are closely fit by curves of the form a^N -|- b. 

Acknowledgments. Special thanks to Lloyd N. Trefethen for pointing out Brut- 
man’s paper [1] and for providing a copy when we had trouble locating the journal. 
Also, we thank a reviewer for pointing out the book m which contains newer results 
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